Spin-current absorption by inhomogeneous spin-orbit coupling 
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We investigate the spin-current absorption induced by an inhomogeneous spin-orbit coupling due 
to impurities in metals. We consider the system with spin currents driven by the electric field or 
the spin accumulation. The resulting diffusive spin currents, including the gradient of the spin-orbit 
coupling strength, indicate the spin-current absorption at the interface, which is exemplified with 
experimentally relevant setups. 
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INTRODUCTION 



Spintronics aims to utilize not only the charge degree 
of freedom but also the spin degree of freedom 1 -. Spin 
currents are an important notion in spintronics from the 
aspect of both basic and applied science. Spin currents 
are generated by a current injection into a ferromagnet 
or spin pumping using the magnetization precession 2 ^ 3 -. 
Other methods to generate spin currents are the spin 
Hall effect 4-6 or the excitation of the spin wave!. In 
the development of spintronic devices, techniques to de- 
tect spin currents efficiently are indispensable as well as 
the spin-current generation. For instance, spin currents 
have been successfully detected via the inverse spin Hall 
effect^—. Due to this effect, spin currents are converted 
into electric currents and therefore are observed as a volt- 
age drop. The electric voltage induced by the spin cur- 
rent has been observed in a lateral junction of a ferro- 
magnetic metal and a nonmagnetic metal with spin-orbit 
coupling (SOC) such as aluminum^ or platinum^ 1 ^. In 
particular, platinum is a typical spin-current detector be- 
cause of its strong SOC due to impurities, which is used 
to absorb a spin current!. 

Kimura et al. demonstrated the spin-current detec- 
tion via the inverse spin Hall effect and suggested that 
the absorption of the spin current occurs from the Cu 
cross into the Pt wire^ -. They explained the spin-current 
absorption using the resistance mismatch for the spin 
current. Namely, they defined the spin resistance as 
R s = X/[aS(l — P 2 )] with the spin diffusion length A, 
the spin polarization P, the conductivity a and the ef- 
fective cross-sectional area S for the spin current, and 
then considered that the spin current flowing in a mate- 
rial is absorbed into the adjacent material with smaller 
spin resistanc e] 10 : 11 . Since platinum is a strong spin-orbit 
coupled material while copper is a weak one, the spin 
resistance of Pt is smaller than that of Cu and hence 
the injected spin current in the Cu cross is partly ab- 
sorbed into the Pt wire. Here, the difference of the SOC 
strength, i.e., the inhomogeneous SOC is the key to the 
absorption of spin current. As shown above, the physics 
of the spin-current absorption induced by an inhomo- 



geneous SOC has been understood only phenomenolog- 
ically, in that previous theories involve the phenomeno- 
logical parameters such as the spin-diffusion length or the 
spin resistanc e) 12 : 13 and microscopic theories of the spin- 
current absorption are missing. The spin Hall effect^ - — 
and the spin-current generation 1 ^ due to an inhomoge- 
neous Rashba SOC have been theoretically examined. 
Also, a spin current generation near an interface due to 
interfacial spin-orbit coupling has been investigated^. 



In the present work, we consider the spin-current ab- 
sorption as follows. In general, the spin-continuity equa- 
tion for spin-orbit coupled systems has a source term: 
+ V • jf = T a . Here, s a , j° and T Q represent the 
spin density, the spin current and the spin torque (or 
source term) with a being a component in spin space, 
respectively. Thus, the divergence of spin currents in 
the steady state is non-zero and therefore leads to the 
generation of spin currents. We refer to a spin-current 
generation induced near an interface between different 
spin-orbit coupled materials or due to an inhomogeneous 
SOC as the spin-current absorption. On the other hand, 
the spin-current absorption in previous phenomenologi- 
cal studies is based on the continuity of a spin current 
at an interface, and therefore absorbed spin currents are 
not generated ones at the interface. 



In this paper, we theoretically examine the spin- 
current generation induced by an inhomogeneous SOC 
due to impurities. We consider the spin-current absorp- 
tion in systems where spin currents are generated by an 
external electric field in Sec. II. In Sec. Ill, we also 
investigate systems where spin currents are induced by 
spin accumulation. We present analytical expressions of 
the spin currents as a response to the gradient of the 
SOC strength. Then, we apply these results to the vicin- 
ity of the interface between metals with different SOC 
strengths, verifying the absorption of spin currents. Sec- 
tion IV is devoted to the discussions. We summarize the 
paper in Sec. V. 
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II. ABSORPTION OF SPIN CURRENT 
DRIVEN BY EXTERNAL FIELD 

A. Model and Formalism 

We consider the ferromagnetic conductor in the pres- 
ence of the spin-orbit scattering due to impurities^— 
and assume that the coupling constant of the spin-orbit 
scattering slowly varies in space. This spatial variation 
of the SOC has not been considered so far. In order 
to address the spin-current absorption, we describe the 
input spin current by the spin-polarized current flowing 
in a ferromagnet under an external electric field. The 
total Hamiltonian is thus composed of the free-electron 
part with the exchange coupling (Hfm), the impurity- 
scattering part ("Himp), the SOC due to impurities with 
its minimal substitution ("Hg + ^so)> an d the interac- 
tion between the vector potential and the electric current 

(^cm) • 

Hfm = E J dr 4{r,t) (~^V 2 - e F ^j c a (r,t), 

(1) 



H imp =J2 I dr U(r) cUr,t)c a (r,t), (2) 

^so = 2- E / dr A soM 

xW(r) • 4(r,t) (V x ow) <v(r,t), (3) 
4 = e ^ / dr Aso(r) 

(7,(7' — ±1 

xVC/(r) • 4(r,t)(A(r,t) x (r aa ,)c a ,(r,t), (4) 

n cm = ^- V [ dr A(r,t)-ci(r,t)V Ccr {r,t). (5) 
2mi J 



a=±l • 



Here, c a (r,t) (cj.(r,i)) represents the annihilation (cre- 
ation) operator of a conduction electron, cr a p repre- 
sents Pauli matrices, e(> 0) is the electric charge and 



(V4) Co 



= 1 and a = — 1 corre- 



spond to the up (f) and down spins (!) , respectively. We 
set epa = e_F + oM with e_p the Fermi energy and M the 
exchange energy. \so( r ) represents the SOC strength 
with spatial variation, averaged over the impurity posi- 
tions. A(r, t) is the vector potential for the external elec- 
tric field and we adopt the fixed gauge as E cm = A. 
U(r) is the short ranged random impurity potential and 
averaging over the impurity positions is carried out as 
(U(r)U(r')) imp = v%n imp 5(r - r'), where u and n imp 
are the strength of the impurity potential and the im- 
purity concentration, respectively. Here, we assume that 
the exchange field is much stronger than the stray field 
and hence we neglect the effect of the stray (magnetic) 
field in our model. Note that, in the Hamiltonian Hg , 
the operator V does not act on Aso( r )- 



The quantum description of the spin current is ob- 
tained by identifying the Heisenberg equation of mo- 
tion for the spin density with the spin-continuity 
equatio n 22 i 23 . The spin current operator in the i- 
direction with a-spin polarization reads 

= IJctfr.tyV^t) - -A lCT Q C t(r,t) C (r,i) 
2m i m 

+Aso E £ajidjU(r)c<(r, i)c(r, t). 
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(6) 

The spin current is thus given by (see also Fig. Q] ) 



1 K 



2m i 



_ e l(k - k '> r (k + fc')itr [a a G kt k>(t)Y 



fc.fc' 



— / ^ Ut E e l{k - k ' +qyr Mq^)t r [*«G k . k ,(t)} 



< 



fc.fc' ,q 

+mJ2 E z l(k ~ k ' +u+pyr ^P \so{u)U(p)ti[G kM (t)]' 



(7) 



Here, G kM (t,t') = ±{T C [c k (t)cl(t')]) denotes the time- 
ordered Green's function (Keldysh Green's function) of 
the total Hamiltonian and G k . k '(t) = lim t /^ t G k . k '(t,t'). 
[• • • ] < means taking the lesser component of the Green's 
functions and tr[- • • ] means taking trace over spin. 

We perturbatively calculate the spin currents induced 
by the inhomogeneity of the SOC strength using the 
Keldysh Green's function formalism. We consider "Hfm + 
"Himp as the non-perturbative part and "Hg j ^so ano - 
"Hcm as perturbative parts. We remark that the non- 
perturbative Green's function contains the self energy by 
the impurity scattering within the first Born approxima- 
tion. Within the linear response with respect to the elec- 
tric field, the second part of the spin currents in Eq.([7]) 
corresponds to an equilibrium spin current, and therefore 
we ignore this contribution. The leading contributions of 
the spin current contain the first order term with respect 
to the SOC strength Aso- Since we are interested in spin 
currents generated by the inhomogeneous SOC strength, 
we focus on the contribution with the spatial derivative 
of the SOC strength. The third part of the spin currents, 
on the other hand, does not involve the spatial deriva- 
tive of the SOC strength. Consequently, the first part of 
Eq.([7]) is relevant to our study. 



B. Local spin current 

We first treat the local spin current, which is driven 
by the external electric field locally. The diagrams of the 
local spin current are shown in Fig. [5] The inhomogene- 
ity of the external electric field is required when we dis- 
cuss the leading effect driven by an inhomogeneous SOC 
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<w>o 

FIG. 1. Diagrams of the spin current for the total Hamilto- 
nian. Bold lines denote the total Green's function. An wavy 
line denotes the vector potential (A(q, Q,)). A broken line and 
double broken line represent the SOC strength (Aso(i*)) and 
the impurity potential (U(p)), respectively. 



where C = -^ri^t v i)- The detailed expressions of 
cti,/3i, and 7, (i = 1,2) are given in Appendix A. For 
a = z, the second terms of Eqs. ©, © and (jTU)) vanish. 
It should be noted that even when M goes to zero, the 
spin currents still have finite contributions and become 
isotropic in spin space. Therefore, the magnetization is 
not indispensable for the spin-current generation itself. 



C. Diffusive spin current 



( a ) \\ S o(u) ( b ) „ ( c ) 




FIG. 2. Diagrams of the leading contribution of the local spin 
current induced by the inhomogeneity of the SOC strength. 
Solid lines denote the free Green's function. Broken lines 
denote the SOC strength (Aso(u)) and double broken lines 
denote the impurity potential (U(p)). Wavy lines represent 
the vector potential (Ai(q,Q,)). 



strength. Namely, the leading contribution of spin cur- 
rents involving the derivative of the SOC strength has the 
form d\sodE cul . The reason is as follows. Spin current 
and electric field are odd under spatial inversion, while 
Aso and other quantities are even with respect to spatial 
inversion. Therefore, the leading terms including spatial 
derivative of the SOC have the forms of d 2 AsoE cm or 
d\sodE cm . We focus on the contributions d\sodE enl by 
assuming <9 2 Aso is negligibly small in this section. Now, 
we present the resulting local spin currents. (For the de- 
tail of the calculation, refer to Appendix A.) Under the 
condition €f<j 3> H/2t (7 (= r/a-) with r a = ?i/27ru§?T,i m pr/ (T 
the spin-dependent relaxation time and v a the density 
of state with spin index a = ±1, we obtain three types 
of the spin currents induced by an inhomogeneous SOC 

0?)^ = C?. a )i 1) + 0'?) i (2) + (7. o, )i 3) : 

(i")f ) = C( ai e a + a 2 (e a x e z )) ■ (VA so (r) x d z E cm (r)), 

(8) 

=C(p 1 (e i xe a ) 

+lh{5 tz e a - 5 ia e z )) ■ VA so (r)(V • E cm (r)), 

(9) 

Crf = C(jie a + l2 (e a x e z )) ■ (VA SO (r) x VE l cm (r)), 

(10) 



(a) \X so {u) 




(b) U so («) 




FIG. 3. Diagrams of the leading contribution of the diffu- 
sive spin current induced by the inhomogeneity of the SOC 
strength, (a) Contributions with the vertex correction to the 
spin-current operator, (b) Contributions with the vertex cor- 
rection to the vector potential, (c) The vertex correction is 
defined as the summation of the ladder diagram with respect 
to the normal impurity potential. 

We calculate the diffusive (or non-local) spin current, 
which can be obtained by including the vertex correc- 
tion to the local spin current. Here, the vertex correc- 
tion means the summation of the ladder diagram with 
respect to the normal impurity potential as shown in 
Fig. [3] (c). Since we focus on the diffusive behavior of 
the spin current induced by the exchange coupling, the 
ladder diagram does not involve the spin-orbit coupled 
impurities. We consider two types of the vertex correc- 
tion, i.e., that for the spin-current operator (Fig. 131(a)) 
and that for the vector potential (Fig. |3J(b)). The vertex 
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correction that cuts across the impurity-averaged line is 
negligible because it is of higher order in r\J£^. For the 
diffusive contribution under the condition tFa^> fi/^T^, 



we obtain two types of the spin currents induced by the 
inhomogeneous SOC (jf)f s ™™ = + (j«)f ) : 



.C 1 ' 



CE a (^e a ■ K^Jr)} +3Me Q x e z ) ■ K^ Jr)}), (a = x,y) 



K[ n) Jr) = I dr> I dr" \ dt' X ^(r -r',r- r",t- t')(V\ so (r>) x 9^ em (r", t')), (n = 1,2) 



(11) 



(12) 



Us)?^ = ^((e, x e«)Ph° + (Si*** ~ <*<«e,)/V) • VA so (r) / dr'^ dt' Xa (r -r',t- t')V • £ cm (r', i')-(13) 



Here. 



X , .,. • (14) 



2vr ^ F CT - iG 



(r 1 ,r 2 ,t)= / —J2 



dn 



2ir jj-< (D a (u + q) 2 - ifi)r ( 



. . r dn^ 



3m 



',(15) 



(16) 



de- 



represent the spin-diffusion propagators. D a 
notes the diffusion coefficient. The detailed expressions 
of a CT , /3i,a-) F a , and Go- (i — 1,2) are given in Appendix 
B. The vertex corrected spin current (jf ) - 1 ' {{jf)^ 2 ') 

corresponds to the local spin current (O'a)i )• 

There does not exist the diffusive-spin current contribu- 

(3) C2') 

tion which corresponds to ■ differs from 

Us )T "* m that only the divergence of the external field 
propagates. 



D. Spin-current absorption 

Let us consider a concrete configuration to show the 
absorption of the spin current by the inhomogeneous 
SOC. We suppose that an external electric field is ap- 
plied to a ferromagnetic metal and therefore a spin- 
polarized current flows in the y-direction as shown in 
Fig. |4] Here, we assume that the applied electric 
field varies smoothly in space from the ferromagnet to 
the attached non-magnetic metal with SOC, which is 
modeled as E cr[L (r) = E y 1+tanh (~ 2 /£g) £y w here £e is 
of the order of the lattice constant. We also assume 
AsoM = ^so@( z )- In the perturbative calculation, we 
have considered a smoothly varying SOC strength. How- 
ever, this sharp spatial variation of the SOC would give 
a qualitatively correct results^ In the present configu- 
ration, a spin-polarized current induced by the external 
electric field flows in the direction parallel to the inter- 
face. Therefore, our setup is different from that of the 
spin-current injection. 



The local spin current in the present configuration is 
generated when Aso an d E cm vary in space, i.e., only 
near the interface between the ferromagnetic metal and 
the spin-orbit coupled metal. Hence, the local spin 
current is not generated in the bulk spin-orbit coupled 
metal. We thus focus on the non-local spin current as 
the spin-current absorption. Since only (VAso( r ') x 
diE ex? {r"))i = E y \ so 5{z'){2i E cos\i 2 {z"/i E ))- 1 5 lz 5 3X 
remains non-zero and V • E evx — 0, the contribution of 
(j£)m (Eq.(ll)) appears in the present configuration. 
From Eq. (|T4"|) . the spin-diffusive propagator in the static 
limit (fi — > 0) reduces to 



-y/B a ID a T a \ri\ 



2.D '(jTij A-u 



ri 



-S(r 2 - n) 



where 



(M 2 + n\f 



M(M 2 



(Af 2 + T)\ ) 3 



(17) 



(18) 



?7 2 7/ + (77 2 (r]l - 3M 2 ) 2 + M 2 {Zril - M 2 ) 2 ) 



x (iMrj+^l - M 2 ) + ia(6r] 2 + M 2 



M 



vl)), 

(19) 

with rj± = . Therefore, we obtain the spin currents 

for each spin component from Eq. (|11[) : 



7T~G 

2£e 



E,.X 



y^SO 



E 



D<jT a 



-Jt 



-i\z\k° 



a„jb„ 



-E y Xso — ° 



26s 



E> a T a 



,-i\z\k° 



A a \JB a 



\AIK 



3 s -- 
Here 



i: 



\y/Ba\ cos(arg(-/Boj) ' 
k° = sin(aig(V^)), 



(20) 

(21) 
(22) 

(23) 
(24) 



which represent the spin diffusion iength and the wave 
number for each spin, respectively. Here, we used the 
assumption b ^> 1° with b the diameter of the interface^. 

We find that spin currents are induced perpendicularly 
to the interface between the ferromagnetic metal and the 
spin-orbit coupled metal, and then decay exponentially 
in an oscillatory fashion as shown in Fig. [4] This ex- 
ponential decay is due to the spin-flip scattering by the 
magnetization in the ferromagnet and also the oscillatory 
behavior originates from the difference between the Fermi 
wave numbers for each exchange spin-split bands. The 
above spin currents are spin-polarized perpendicularly to 
the spin polarization of the input spin-polarized current 
(namely, z component). This reflects the violation of the 
conservation law of spin near the interface, and in other 
words, spin torque generated at the interface produces 
the present spin-current absorption. 

It is also found that the magnitudes of the resulting 
spin currents in Eqs. ([20]) and (j2"Tj) at z — are, respec- 
tively, proportional to rj/ and?/_ry + in the dirty limit 

(t ct — > (a =J, J,))) and r]-r]V 2 and rj^ 2 in the clean limit 
(r a — » oo (<r =t,|)). Namely, the magnitude of the ab- 
sorbed spin currents gets reduced in clean ferromagnets. 
This is due to the competition between the two effects of 
impurity scattering. In general, a SOC due to impurities 
in clean ferromagnets is weak due to low impurity con- 
centration, leading to a small spin current. On the other 
hand, low impurity concentration leads to a large spin 
current since impurity scattering is suppressed. When 
the former effect overcomes the latter, spin currents are 
expected to be suppressed in clean ferromagnets. Finally, 
we note that for junctions composed of a weak spin-orbit 
coupled ferromagnet and a strong spin-orbit coupled non- 
magnet (namely, Aso( r ) = A so 6>(z) + ^so^(~ z ) wrtn 
Ag Q > A| ), spin current absorption predicted in this 
paper also occurs. 



Non -magnet A 50 > 



III. ABSORPTION OF SPIN CURRENT 
DRIVEN BY SPIN ACCUMULATION 




Ferromagnet A so = 



FIG. 4. (Color online) Schematic illustration of the spin- 
current absorption by the input spin-polarized current. The 
external electric field E cm is applied to a ferromagnetic metal 
and therefore a spin-polarized current j* flows in the direc- 
tion of y axis. The non-magnetic metal attached to the fer- 
romagnetic metal has a strong SOC Aso due to impurities. 
The diffusive spin currents jf' y are induced perpendicularly 
to the interface and then decay exponentially in an oscillatory 
fashion. 



"Ho j W-a.cc, Himp and Hg . Here, 

U acc = ^2 / dr MM) cl(r,t)c a (r,t) 

= E J dr cl(r : t)(fl(r,t)S aa/ + n s (r,t)a* a ,)c a >(r,t). 

er,<x'=±l 

(26) 



(25) 



Here, we have defined the spin accumulation as the differ- 
ence between spin-up and spin-down chemical potentials, 
namely, jit s (r) := PT ^ ; 12 i 13 This spatial distribution of 
spin accumulation leads to a diffusive spin current. In 
the following, we will also assume the chemical potential 
p(r,t) to be constant. 



B. Spin-current 



A. Model 



In this section, we examine the spin-current absorp- 
tion originating from the variation of the SOC strength 
in systems with spin accumulation. These systems are 
modeled by the conducting electrons in the presence of 
the gradient of the spin-dependent chemical potential and 
that of the SOC strength. Our Hamiltonian consists of 



We will calculate the spin current corresponding to 
the diagrams in Fig. [3] in a similar manner to the pre- 
vious section. Here, we have replaced the vector po- 
tential A(q, fi) with the spin accumulation fi s (q,Q) in 
the diagrams of Fig. [31 We consider Ho as the non- 
perturbative part and T-L aC c, "Wimp and Hg as pertur- 
bative parts. In this section, we adopt the self energy 
including both the normal and the spin-orbit coupled 
impurity potentials'.^. We also treat only the vertex 
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correction for the spin-current operator—. This correc- 
tion provides a finite decay length of the spin current 
determined by the SOC strength. After taking the lesser 
component of the Green's functions and averaging over 
impurity positions, the spin current shown in Fig. [2] (a) 
and (b) reads 



UfU = i-h 

m 



2?r' 



VLe tnt 



E E 



i,j,k h,h' ,u.q 

xe^ u+ ^ r X so (u)^(q,n)(U(p)U(~p)) imp 
Jk 



xeijfcR(Xc«jfcm)tr [<r a (7 k <7 z ] 



(27) 



For the expression of Xaijkmi refer to Appendix C. In 
the same manner as in the previous section, we expand 
the coefficient Xaijkm with respect to u and q. Ac- 
cording to the inversion operation r —> — r, it is found 
that only the contributions with an odd order of u 
and q remain. Now, we perform the vertex correction 



Xaijkm^ 



Xaijk Z^a,b,c,d,e a da L ad^cb^be 



with Xaijk being the vertex-corrected coefficient and 
F a d,cb being the vertex correction corresponding to the 
ladder diagram including the normal and the spin or- 
bit coupled impurities^. The contributions from the 
vertex-corrected coefficient with the first order of u or 
q vanish, and hence the leading contribution in the dif- 
fusion regime is of the third order of u and q: Xaijk — 

-(2,1) -(1,2) 
Xaijk ' Xaijk 



where the superscript (i,j) denotes 



the order of u and q. The contributions from xfiik ano - 



where £(i,2) 



aijk 

Xaijk vanish exactly. 

(1 2) 

From the contribution of Xaijkm' we obtain the expres- 
sion of the spin current: 

(j})£ 2) = C(i,2) J dr> J dr" J dt' D(r -r',r- r", t - t') 

x [d a X so (r') d z (d mf i s (r",t')) 

-d z Xso(r') d a (d mt M s (r",t'))}, (28) 
I 



with t = h/2i 



3m ,viuii ' — "/ a O"imp 

relaxation time and v being the density of state 



v being the 



D(r u r 2 ,t) := / — £\ 



^i£lt-\-iu-ri -\- iq-r<2 



n(i -*)(! + 3k) 
Dr((u + qf + e) + Mr V ' 

is the spin-diffusive propagator with k = Ag /c|,/3 (fcp- 
is the Fermi wave number), D = 2ci?r/3m and £ = 
d 4k/ Dt. As seen from Eq.(29), £ determines the de- 
cay length of the spin current. 

As for the other contribution Xaijkm ■ 
spin current of the form: 



we obtain the 



(2,1) 



C(2,l) 



x [d af x a (r",t') d z (d m X so (r')) 



dr' / dr" / dt' D{r -r',r- r", t - t') 



-d z fi s (r",t') da(d m X so (r'))} , (30) 

where C( 2 ,i) = 60 % h m TV (nr) • Here ' we have focused 
on the diffusive spin currents, corresponding to contri- 
butions with the vertex correction. It is found that if 
the scales of spatial variations of Aso( r ) an d ^( s (r,t) 
are nearly the same, the magnitude of the spin current 

~ (2 1 ) ~ ( 1 2^ 

{jf)m is larger than that of (j")m by the factor of 
cft/H. 



C. Spin-current absorption 

Let us consider a typical configuration in order to in- 
vestigate the spin-current absorption. We assume that 
two non-magnetic metals with different SOC strengths 
are connected at the z = plane as shown in Fig. [S] In 
addition, there exists a gradient of the spin accumulation 
along the y direction in one of the metals. In the vicinity 
of the z = plane, we assume that Aso( r ) = Aso(z) and 
[i s (r,t) = /j, s (y,z,t), leading to 



(if), = o, 

(jf) z = J dr' J dr" J dt' D(r -r',r- r", t - t') 

x [-C(i,2)d z "Aso(2") d zl (d y/f i s (y',z',t')) + C {2A )d y/ fis(y\z',t') d 2 z „X so (z")] 
(It)* = 0. 



(31) 

(32) 
(33) 



The first term of the right hand side in Eq. (f32|) rep- 
resents the interplay between the variation of the SOC 
strength and that of the input spin current near the in- 
terface. Note that spin current is induced by the gradient 
of the spin accumulation without SOCJ^ On the other 
hand, the second term comes from the variation of the 
SOC strength and uniform input spin current. Remark 
that a spin current polarized along the y-axis is generated 



as a response to that polarized along the z-axis. Similar 
to the absorption of spin currents driven by an external 
electric field, a spin current absorbed into the top metal 
is non-local, and diffusive spin current decays due to the 
spin-orbit coupled scattering as seen from Eq.(29). 
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Non-magnetic metal with SOC 




Non-magnetic metal without SOC 



spin currents in our theory are of different origin from 
those in the theory based on the spin-resistance mis- 
match. Hence, we cannot compare our results with pre- 
vious results directly. Both studies show the spin-current 
absorption with different mechanisms. 

In asymmetric structures, conduction electrons at its 
interface in general feel a Rashba SOC 2 ^, leading to the 
spin polarization under an external electric field or cur- 
rent injection 26 . In the configurations of Figs. [4] and 
03 there would exist a Rashba SOC at their interfaces. 
The current-induced spin polarization due to the Rashba 
SOC may affect the absorption of spin currents predicted 
in our models. 



SUMMARY 



FIG. 5. (Color online) Schematic illustration of the absorp- 
tion of spin currents driven by the spin accumulation. The 
arrows in the bottom metal denote input spin currents (j* j n ) B 
by the spin accumulation. In the top metal, the left arrows 
are normal diffusion of the input spin currents (js- ln )z, which 
is irrelevant to the variation of a SOC, whereas the middle 
and the right arrows are diffusion of spin currents (j'" abs )z 
induced in the vicinity of the interface, relevant to the differ- 
ence of SOCs. 



IV. DISCUSSION 



We have investigated the generation of spin currents by 
an inhomogeneous SOC due to impurities, which have 
been applied to the interface system to show the spin- 
current absorption. Using the Keldysh Green's function 
formalism, we have presented analytical expressions of 
the spin currents with the gradient of the SOC strength 
for two systems: the system with field-driven spin cur- 
rents and one with spin-accumulation-driven spin cur- 
rents. The resulting spin current indicates the absorp- 
tion of the spin current at the interface between materials 
with different SOC strengths. 

In the present study, we assumed a homogeneous mag- 
netization. The extension of our model to an inhomoge- 
neous magnetization is also an interesting future work. 



In the present study, we have focused on diffusion of 
spin currents including the gradient of the SOC strength. 
However, there also exists a diffusion of spin currents ir- 
relevant to the difference of SOC strength as shown in 
Fig. [5] We can extract the former contribution by com- 
paring spin currents in the triple lateral spin valves us- 
ing middle junctions with different SOCii. In fact, the 
absorption of non-local spin currents has been success- 
fully demonstrated in the spin-valve measurement con- 
sisting of double Py/Cu junctions and middle Au/Cu 
junctio n 10 ' 11 . 

Since the argument based on the mismatch of the spin 
resistance by Kimura et alr^ assumes the continuity of 
a spin current, the absorbed spin currents are not those 
generated at an interface. Consequently, the absorbed 
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Appendix A: Calculation of local spin current 



We perturbatively calculate the spin currents induced by the inhomogeneity of the SOC strength using the Keldysh 
Green's function formalism. We first treat the local spin current, which is locally driven by the external electric 
field. The leading contribution of the spin current involves the first order with respect to the SOC strength, which 
is diagrammatically shown in Fig. [5] The local spin currents with a-component spin polarization flowing in the m 
direction, i.e., (j?) l ° cal = + thus read 



?(.■»+<!>*• \ so ( u )(U(p)U(- P )) imp Mq,n)e ijk 2$t( X aijkl m ), (Al) 



i,j,k,l k,k' .u.q 

{3t)t = i^e ^fie int £ E e i ^-'-Aso(«)(^(p)^(-p))i m p^(9^)^ fe 25R( Xmfcm ), (A2) 

i,j,k k,k' : u.q 

where the superscripts "sj" and "sk" denote the side-jump and the skew-scattering contributions, respectively^ and 
each coefficients are given by 

Xaijklm =+(k-u~ k')i(k + k')j(2k -2u- q)i(2k -u- q) m tr [a a g^a k g§,g§_ u g^_ u _ q ] 
+(k' - k)i(k -u+ k%(2k -2u- q)i(2k -u- q) m ti [a a g£g%a k gt u 9k- u - q ] 
+(k-u- k%(k + fc') J (2fc / - q)t(2k -u- q) m tv [a a g£o- k g*g£,_ q gt u - q ] , (A3) 
Xaikm = (k-k' -U- q) l (2k -u- q) m tr [a a g^a k 'g£>9k- u -q\ ■ (A4) 



Here, g k = 9fc L=o> an d 9k u an< ^ 9k u denote the non-perturbative retarded and advanced Green's func- 
tion, respectively, which are 2x2 matrices in spin space. We also use the Langreth's method to obtain 
' ' •fl'fei,u;5fe 4+ i,u>-n • • •3fc„,c J -f2] < ^ -^f( u )9ki, u ---9k i , u 9k i+1 ,oj---9k n , u with the Fermi distribution 

function^. Here, we set temperature to zero, i.e., f'(cj) — —6 ((J). 

Since we are interested in spin currents produced by the inhomogeneous SOC strength, we focus on the contribution 
which involves the spatial derivative of the SOC strength. After some algebraic calculations, we obtain the above 
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coefficients represented as 
2 

Xaijklm ^k 2 (S ]lUl q m + 5 ]m u t qi) (tr [a Q g*a k g%g%g£] + tr [a a g% g%(J k g% g£ ] ) 
4 ft 2 

+ T 7-fc 4 (* 3 !« ! 9m + <W*i?i + *m«i«f) (tr cr fc $#s£ (ffjf ) 2 ] + tr [a Q 5 ^a fc 5 «( 5 ^) 2 ]) 



2 

+ 3 (k' 2 SjlUiq m + k 2 S ]m u t qi) tr [a a g k a k g k ,g k ,g k 

o fc2 ^ fc2 



9 m 



Xmfem ~ (w l9m + qi u m ) ( tr [<J a g%a«g%g%\ + —k 2 tr {<T a g£<J k g&(g£)^ 

m 



(A5) 
(A6) 



Here, we assumed the rotational symmetry of the system, i.e., J2k = J2k T^J anc ^ Sfc hkjhk m — J^k fs 



$u&jm + Si m Sji). We remark that if we do not consider the inhomogeneity of the external electric field, the present 
contributions vanish. The inhomogeneity of the external electric field is required to obtain the contributions from the 
inhomogeneous SOC strength. 

Next, we take traces over spin space in Eqs. (|A5I) and (|A6|) by using the formula ti[a a Aa k B] — X^=±i(^" fc — 
SazSkz — ivtazk)A a B s + J2a=±i ^az^kzA a B a with A, B Green's functions, and then Eqs. (|A1[) and (IA2|) reduce to 

ia,<Z i,j,k,l a 



-5jiUiq m I (loi;00 + ^10;Oo) ( Qll-,01 + k^H;02 ) + QoijOl/lOjOX + oQoiiOlQ 



'01:00 + ^10;Oo) 



15 

eft 2 
~2^ 



« (ft 

filmUiqj (Iqi-oo + -fl0;0o)'^ll;02 



Qll;01 + g^ll;02 



7oXi01^10;01 + 3 < 9oi;02 ( 3lO;01 



it,q i,j,k,l 



5jiUiq m oo ^Q^oao + g^20;20 J + ' 



?10;loAo;10 + 3 QlO;loQlO;20 



--^SjmUiqi ( 2/f 0;00 ( Q21 



111:00 I l v'2(i:IO + g^20;20 



^10;loQlO;10 



o Ql0;20*3l0;10 



16 



+ -^g^im^i^jAojOO^O^O 



(A7) 



i,j,k u.q 



^2(Sak ~ S az S kz )?ft + <J€ ak Q) 



xu iQm,IoO;Ol(IlO;01 + ^QlO-m) + S az 6kzUiq m ^tIoO;10 (Ao;10 + 2Q 10;2 o) 

(X 

Here, we define integrals of the Green's functions appearing in the above expressions as 

Iab;cd = y~!(gfc.g) a (fffc i g) b (gfc,g) C (gfc.g) d ; 
k 

Ql b , cd ^Y, ( u(9l.) a {9U\ 9 L) c {gi s ) d , 

k 

Sab:cd — ^ e k(9k,a) a (9k,a) b (9k,a) C (9k,a) d J 



(A8) 

(A9) 
(A10) 
(All) 



where g% a = [tFo-tk+iVo] 1 with rj a = h/2T a and g£ a = {g k ^)* with e k = We sum up Eqs. (|A"7|) and and 

consider the dominant contribution from the integrals of the Green's functions under the condition e^o- 3> h/2r a . By 
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carrying out the integrals of the Green's functions in Eqs.(A9-All) and transforming to the real-space representation, 
we obtain the local spin current of the form 



■net? 
2m 



(vf + vi) [(aie a + a 2 (e a x e z )) ■ (VAso(r) x d m E cm (r)) 

+(0x(e m x e a ) + ^2{S mz e a - 5 ma e z )) ■ VA S o(r)(V • E cm (r)) 
+(7ie Q +72(e Q X e z )) • (VA so (r) x VE^(r))} . 



Here, each coefficients are dimensionless and given by 

^ = 7i + Si, 

(- 16 / e F M 2 r,+ _ 1_1 rj+ y, \ , _ s 



5? = — 



(a = zj 

16 / e F M{M 2 -rf + ) 1 1 



M 



9 ^ (M 2 +r? 2 ) 2 4 1^ + n M 2 + J7 2 



£2 



7i 



A/ 



9 VM 2 +r )+ 

_^ CT 

9 ?7 

»7+ v- < /•■ 



9 M 2 + rf\ 



_2_ 

15 



-A^, (a = z) 

15 1)+ ' v / 



(a = x,y) 



72 = -77 



M 2 (?? 2 -jj 2 ) 



(M 2 + r; 2 )(Af 2 +T/ 2 ) 



V+(M 2 — r}\) 



e Fa t e F M 



(M 2 + 7]\) 2 ^ V* M 2 + 7 1 \ 



(A12) 

(A13) 
(A14) 

(A15) 

(A16) 

(A17) 
(A18) 

(A19) 

(A20) 



and rj± 



Appendix B: Calculation of diffusive spin current 

We calculate the diffusive (or non-local) spin current, which is represented by the vertex correction to the local spin 
current. This can be performed by replacing Xai 3 ki m and Xaikm in Eqs. ([AT) and (|A2]) with Xaijklm = ^u, q n -u, q Xaijkl + 
r q n gXSjfem and Xaikm = ^q^gXaik, respectively, where 



K {u) = ^pvy {2k _ ^Xmmw. (Bi) 




H*r(«) = ( =^ Y^g^gt-^-n, ) ■ (B2) 

n=0 

(B3) 



X%u = +(*; - u - fc')i(fc + fc'),(2fc - 2u - 9 ),tr [<T a g%<T k g%g£_ u gtu- 

+(k' -k)i(k-u + k') 3 (2k - 2u - ? ),tr [<T a g£g%<T k g£- u gLu-q\ 

- « - fc')i(fc + fc'), (2fc' - «),tr [a a g£a k g%g&- q gL u - q ] , (B4) 
= +(fe - « - fc')i(fc + (2fc - « - <?) m tr [* a gZ<T k g%gL u gt 
+(k' - k)i(k -u + k') 3 {2k -u- q) m tr [<T a g%g%<r k g£- u gt 
+(k-u- k%(k + k%(2k -u- q) m tr [<T a g^ k g^_ q gt u - q ] , (B5) 
Xmfc = (fc - k' - u - q)itr [<J a 9%<T k g£,gt u - q ] ■ (B6) 



-q\ 



-u—qj 

-u—q\ 
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Expanding with respect to u or q, we obtain the leading contributions as follows 



(i) 



(2) 



Xaik 



-\kH m (tr [<j»g* o k g%g%g£} + tr [^S&S^M ) - -fc^tr [a a g%* k g% g& g£] 
a n 2 

-—kVs^u + g) it r [* a g£* k gZg£(g£) 2 ] 



-lk% mUi (tr [a a gZa k g%g%g£] + tr [a a 9 M^ 9 M]) ~ f^i^tr [^a fe 9 £«^] 
4 fi 2 

: -(« + g), (tr [^a^j*] + i-fc 2 tr [^^(gj*) 2 ] 
\ 3 m 



(B7) 

(B8) 
(B9) 



Next, we take traces over spin space in Eqs. (|B7[) . 
then the diffusive spin currents are reduced to 



and (|B9[) in a similar manner to the local spin current, and 



u.q i,j,k,l a 



'ak - SazSkzW + cre ak ?s) 



"^SjlU%q m H a ,a(u + q)u ni mp I 10 . 01 ( (1q1;00 + ^?0;Oo)Qll;01 + Qoi;01^10;01 + g Qoi;01 ( 3lO;02 
+ ^jm u iQl n <JjCr (g)MQni m p/J r 0;10 ((/o 1;0 o + -^10;Oo)'9ll;Ol) + A)l;0lQl0;01 



3 

2m 



w,q i,i,k,l a 



— SjiUiq m IT^o- (it + g)M ni mp /j T . 1 Q I 2/f . 00 Q20;10 + Ql0;10^10;10 + gQlOjloQlO; 

2 

+ g<5j m Ujg; n <TjO -(g)it ni m p/ 1 Q ;10 (2/ 10 .Q <520:io + ^loaoQiOjio) 



^(5 Qfc 5R + ae ak ^s)(uiq m + q l u m ) U a ^(u + q)uln imp I° . 



01 \ A)0;OlAo;01 + g^00;01 ( ?10;02 



' y ] 8 az $kz$t(Uiqm + QiU m ) Hcr :<T (u + gjuo^imp -ffo;lQ ( ^00:10^10:10 + g ^00;10 ( ' ! 1 1 : _» ! 1 



(BIO) 



(Bll) 



Here 



n CTl(7 (tt) 



cr 2 cr ^ 1 2 2 

1 — ( Wo^impAo^l "I" M 7l i m P-^10:02 ^^2m 3" U nim pQl0;03 

ft 2 1 

1 — ( Mgnimp/jQ.jo + Ml UQni mp I 10 . 2 Q + ^2^ 3^ Uoriim P^ 10 ' 30 



(B12) 
(B13) 



By carrying out the above integrals of the Green's functions and transforming to the real-space representation using 
u\ so (u) = -ij dr' V r ,\ S o(r')e~ lr '- u and Qq m A(q,fl) = - J dt' J dr" d' t ' n A(r" , t')e~ ir " - q e- int \ we obtain the 
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diffusive spin current of the form 



{Sax + Say) ^2{e a Ua a + (e Q x e 2 )3a ff ) ■ J dr' J dr" J dt' 

(T 

/ h E (VAst>(r,) x 

+5 Q2 Y^a a e z - J dr' J dr" J dt' 

r fin _ „iO(t— t')+*u-(r— r')+»g-(i — r") 

/ IftMtf-M),. (VAso(r-)xa„ E .„ (r ", t -)) 



/• /• /■ e tfl(i-t')+<g-(r-V) 

^(/3i, CT ( ei xe Q )+/3 2i(T ((5 l2 e Q -5. iQ e 2 ))-VAso(r-) y dr'y dt' J — ^ — g2 _ — - V-E cm (r',t') 

(B14) 



where 



2 M„ + ((l+a)r£-M 2 ) 2 2 ^((l+q)M 2 -^) 2 , _ , 

9 (Af 2 +r; 2 ) 3 e -FV T fg (A/2+^)3 e i?o-! l^CK — J,, 



-I<L£v ( a = z) 

»*2 , 2 2 

2 
3 



(B15) 



ft (M 2 +r, 2 )(M 2 +r f 2 ) + 71 ft "fl r,,, > V x ~ ^ : U ) (B16) 

^^nn^. (« = *) 

P2 ' CT "3 ft (Af 2 + 77 2 )(M 2 + 77 2 ) ' 1£>± ' J 

^ _ M 2 2qMr^ 2 t? 2 ^ 2 (ry 2 - 3M 2 ) 2e Fff T CT 2 

^ = + (M 2 + ^)2 "^ + ( M 2 + ,2 )3 3^T ( " + ^ T " (B18) 

Gct = MHi + (M 2 + , 2 ) 2 " Tg + (M 2 +r; 2 )3 3^ ( " + q) T " { ™ ] 

D a = ^. (B20) 
dm 

Appendix C: Calculation of spin current driven by spin accumulation 

The spin current shown in Fig. [5] (a) and (b) reads 
°'" )m = / £ ^ ? 51 e l ^-'-AsoN(C/(p)C/(-p))i m pM s (q, tt)e iifc tr 2R( Xaiifcm ) 

k,k' ,u,q 

fa? f d 

= " 4 ^ ftu n imp / 2^X! 51 el(M+<j) ' ? " A S0(M)M ;5 (9,^)('5 l2 ^a -^^^^(Xmifcm), (CI) 

Xaijkm = {k-u- k')i{k + k') :j {2k -u- q)m9k9k'9k-u9k-u- q + i k ' ~ fc M fc - u + k ')j( 2k ~u- q)m9k9k'9k-u9k-u- g 
+{k-u- k'Uk + k'),{2k -u- q)m9%9k-9k>- q 9k-u- q - (02) 
As for the above coefficient Xaijkm, the contribution involving two q's and one u is calculated as follows; 



XaifL = ~oTT- [2k 2 u iq , nqj g%{9k?{9£? + k^q^g^g^g^ig^) 2 + k n u lWm g*gl{g£,f g£\ 



2 ft^ 

3 2m 

l (£) 2 flSfo&) 2 «?Cfl£) 9 - (03) 
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The resultant spin current reduces to 

*2 



Ci. tt )& ,2) =-4^^n lmp J ^ ne im Y / ^ u+qyr Xso(u)fi s (q,n)(u a q m q z ~u z q m q a ) ^ [l? . 00 Q% 0;20 + /f 0; i Qio ; 2o] 



(C4) 



Now we perform the vertex correction Xaijkm^ [a a a k a z ] -> Xoijfc S a ,6,c,d,e CT da r ad,cfc°6e 'f c witn r ad, c fc being the 
vertex correction corresponding to the ladder diagram including the normal impurity and the SOC due to impurities^, 
and 

Xaijk =(k-u- k')i(k + k')jg%g§g%_ u g£_ u _ q + (k' - k)i{k -u+ k')jg%g%g%_ u g£_ u _ q 

+(k -u- k%{k + k') jg Xg*g£_ q gt u - q . (C5) 

According to Ref. the vertex correction is calculated as T a d. c b{q, fi) = Tc(q, Cl)6 a dS c b + ^s{q, H) J2i a ad a cb w ^ n 

1 + 3k 



r c (g,n) = 
r s (9,n) = 



2(ZV +^)t' 

(1-k)(1 + 3k) 
2(4k + (Dg 2 + in)r) : 



(C6) 
(C7) 



leading to 



^adMe^c = r c tr[a Q i]tr[ia V] + £ r s tr[<r Q £ r , ]tr[a , aV] 



,,b.c.d,< 



4ie Q fe 2 rs. 



(C8) 



Here, r = fi,/27TUQni mp ^, k = Ag fc^/3 (fop is the Fermi wave number), D = 2ei?r/3m, and £ = d '4k /Dr. 1 is the 
identity matrix. Therefore, we obtain the spin current including the vertex correction of the form 

4-kT?tv 



(j's )m 



2) 



3m 



dr' / dr" / dt' 



2^ 



E- 



,in(t-t')+»tt-(r-V)+ig-(r-r") ^/ r (l ~ + 3«) 

D((u + q) 2 + e) + iSl 



x (d a X so (r') d z (d mf i s (r",t'))-d z X so (r') d a (d mf i s (r", *'))) , 

with ^ being the density of state. We remark that the spin-diffusion propagator is given by 

frfJ]v^e ,!!t+ra ' ri+ *« T2 Sl/T(l-K)(l + 3K) 
D( ri ,r 2 ,t):= — ^- 



(C9) 



£((M + q) 2 +£ 2 ) + zft 



= — 5()"2 — ri)^ 



2 N 2 i 

4L>£ 2 2f 



exp 



2 , , KI 2X1 



4L>t 



(CIO) 



On the other hand, the contribution involving two u's and one q is given as follows; 



Xaijlm = -^^ k2u * u j<irnlgkgk>gk) 2 g£ ~ ^ (^) k 4 (uiUjq m + UiU o q S jm + UiU m qj)gk9k>(9k) 2 (9k) 2 
4 h 2 l 

-g —tfuiiumqj + q m Uj)g%g%g% (g£ ) 2 - -k 2 { Ul u mqj + u lU:j q m )gj}g§g£(g£) 2 
~^k' 2 u l u m q ] g^g§{g^) 2 gt - ^k 2 k' 2 u t u m q 3 g§ g% (g£) 2 (g£) 2 . 



(Gil) 



The corresponding spin current reduces to 

U:)t 1} = -4^-hu 2 n inip f ^ m nt J2e l{u+q) - r Xso(u)^(q,n)Y, e ^ e ^ 



u.q 



4 <r 4 CT .... 

g Ao;00 ( 920;20 — g Ao;10*9lO;20 g~(QlO;20 



a \2 
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TtJ OCT 

J 10;00°30;20 



^i^oqo^jm^- 



1G 
15 



TCT Off 

J 10;00 J 30;20 



(C12) 



By performing the vertex correction, we obtain the final expression of the spin current 



152itH 2 tv (2e F T 
45m \ H 



iO(t-t')+w-(r-r')+ig-(r-r") ^/ T (l ~ + 3«Q 



u.q 

x (a a M s (r',t') a,(ft n Aso(r"))-fl*M-(»'',*') 9 a (5 m A so (r"))) . 



D((u + q) 2 + e) + ifi 

(C13) 



